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We study the zero-temperature phase diagram of a gas of bosonic 87 Rb atoms in two-color superlattice po- 
tentials starting directly from the experimental parameters, such as wavelengths and intensities of the two lasers 
generating the superlattice. In a first step, we map the experimental setup to a Bose-Hubbard Hamiltonian with 
site-dependent parameters through explicit band- structure calculations. In the second step, we solve the many- 
body problem using the density-matrix renormalization group (DMRG) approach and compute observables 
such as energy gap, condensate fraction, maximum number fluctuations and visibility of interference fringes. 
We study the phase diagram as function of the laser intensities S2 and s± as control parameters and show that all 
relevant quantum phases, i.e. superfluid, Mott-insulator, and quasi Bose-glass phase, and the transitions between 
them can be investigated through a variation of these intensities alone. 



Ultracold atomic gases in optical lattice potentials offer 
unique possibilities for studying fundamental quantum phe- 
nomena, such asphase transitions to exotic quantum phases 
fl, & & II H U 01- Experimentally, the key advantages 
of these systems are the unparalleled in-situ control over all 
relevant parameters and the flexible tools to probe the quan- 
tum state. This has lead to a series of experimental studies, 
which recently also focus on non-homogeneous lattices and 
two-color superlattices 0, H, S 0] • Theoretically, these sys- 
tems are well described by Hubbard-type Hamiltonians and a 
whole range of approaches from mean-field many-body meth- 
ods to exact diagonalization and density-matrix renormaliza- 
tiongroup techniques is being used to study the phase diagram 
tfH la EflUll [12, IS [3 Theoretical studies typically 
adopt the generic parameters of the Bose-Hubbard Hamil- 
tonian as control parameters, which entails several simplifi- 
cations as compared to the experimental situation. A zero- 
temperature phase diagram for a Bose gas in a two-color 
superlattice in terms of the generic Hubbard parameters in 
shown in Fig. Uta). 

In this work we want to establish the direct link between 
the experimental setup as, e.g., discussed in Refs. JUS] and 
the phase diagram by using the simplest experimental vari- 
ables, the intensities S2 and si of the two lasers generating the 
superlattice, as control parameters. To this end, we first de- 
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FIG. 1: (color online) Contour plots of the energy gap computed 
in DMRG for a commensurate superlattice with / = N = 30 as 
function of the generic Hubbard parameters (a) and a function of the 
experimental laser intensities (b) (see text). The labels mark the do- 
mains of the superfluid (SF) phase, the homogeneous Mott-insulator 
(MI) phase, and the quasi Bose-glass (BG) phase. 



termine site-dependent Hubbard parameters from the experi- 
mental superlattice setup using explicit band- structure calcu- 
lations. We then solve the many-body problem using state-of- 
the-art density-matrix renormalization group techniques and 
map out the phase diagram as function of the laser intensities. 
An example for such a phase diagram is shown in Fig. [TJb) 
and will be discussed in detail in the following. 

Bose-Hubbard Hamiltonian & Band Structure. Strongly 
correlated ultracold bosonic atoms in a sufficiently deep one- 
dimensional optical lattice potential with I sites are well de- 
scribed by the single-band Bose-Hubbard Hamiltonian fl6ll 



H 
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Here a J (a f ) are creation (annihilation) operators for a boson 
in the lowest Wannier state localized at site i and n i = aja^ 
are the corresponding occupation number operators. The three 
terms of the Hamiltonian describe the tunneling between ad- 
jacent sites, the on-site two-body interaction and the on-site 
potential, respectively. The site-dependent Hubbard parame- 
ters J^i+i, Ui, and €i, which control the strength of the indi- 
vidual terms, contain all information on the underlying lattice 
potential and the interaction between the atoms. 

In most previous applications, the Hubbard parameters 
were used directly as the control parameters spanning the 
phase diagram. Superlattice systems are typical approximated 
by using site-independent tunneling and interaction matrix el- 
ements, i.e. Jj,i+i = J and Ui = U, and some ansatz for 
the site-dependent on-site energies ei (THEHl. I n this study 
we start directly from the experimental setup, i.e., from a one- 
dimensional lattice potential V(x) in real-space depending on 
the wavelengths Ai and A2 and the intensities s\ and S2 of the 
two standing waves through V(x) = s\E ri sin 2 ( + </>) + 

S2E r2 sin 2 (f^-^), where E ri = 2 m\ 1 * s me reco ^ ener gy of 
the atoms with mass m and (j) a relative phase shift between 
the two standing waves. For the sake of simplicity we do not 
account for a longitudinal trapping potential. We study two 
different superlattice topologies: A commensurate superlat- 
tice with A 2 = 800 nm, Ai 1000 nm (<j) 7r/4) and an 
incommensurate superlattice with A2 = 830 nm, Ai = 1076 
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nm ((f) = 7r/3). These values are motivated by the experimen- 
tal setup discussed in Ref. J3l and we also assume the strong 
primary laser to be defined by A 2 and 82. 

In order to extract the Hubbard parameters for the lattice 
potential V(x) we first have to determine the localized Wan- 
nier functions via a single-particle band- structure calculation. 
For a homogeneous lattice (si = 0) with / sites we can eas- 
ily obtain a numerical solution for the Bloch functions of the 
lowest band, ^ (x), characterized by a quasimomentum index 
fc = 0, 1, — 1. The corresponding Wannier functions are 
then determined by a Fourier transformation with respect to 
the quasimomentum, Wi(x) = -^j Ylk=o e~ Y ^ kl e 1(pk ipk(%)- 
The arbitrary phase cfk of the individual Bloch functions is de- 
termined by requiring the Wannier functions to be localized. 
From these Wannier functions the Hubbard parameters can be 
computed through matrix elements of the different terms of 
the real- space Hamiltonian 

f ( Ti 2 d 2 \ 

-Ji,i+i = dx w*(x) ( ~ + V ^ ) 

6i = j dx W*{X) (-^^2 + V ( X )) W ^ X ) (2) 

Ui = 2hu± a s J dx \wi{x)\ A . 



The interaction is described by a three-dimensional s-wave 
contact interaction with a strength given by the s-wave scatter- 
ing length a s . For the transverse directions the wave function 
is approximated by a Gaussian of a width corresponding to the 
frequency co± of the transverse confinement. Following Ref. 

we consider a gas of 87 Rb atoms with s-wave scattering 
length a s = 109a Bo hr and we assume a transverse trapping 
frequency 00 ± = 30E r2 /h to ensure the validity of the ID 
description. 

In the case of non-homogeneous lattice potentials we de- 
scribe each site individually with a Wannier function extracted 
from a homogeneous lattice with the same local depths. Based 
on this set of nontrivial site-dependent Wannier functions the 
Hubbard parameters are then determined from Eqs. ©. This 
scheme is sufficiently accurate for the parameter regime under 
consideration. Even simpler approximations, e.g. a lowest- 
order perturbative inclusion of the lattice inhomogeneity by 
computing the Hubbard parameters © with the full potential 
V(x) and the Wannier functions from a homogeneous lattice 
of the same average depth, yield similar results. 

Examples for the site-dependent Hubbard parameters for 
the commensurate as well as the incommensurate superlattice 
with S2 = 10 and s\ = 1 are depicted in Fig. [2 Note that we 
always subtract a global constant from the Hamiltonian such 
that e m i n = 0. The dominant effect of the superlattice struc- 
ture is the spatial modulation of the on-site energies which 
is in-line with the approximation to introduce the superlattice 
through €i only. However, also the tunneling matrix element 
J^i+i, which essentially depends on the barrier height be- 
tween the sites i and i + 1, shows a sizable variation of ±20% 
around the average value J. The interaction strength Ui shows 
only a weak site-dependence which is induced solely through 
the site-dependence of the Wannier functions. A comparison 
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FIG. 2: Site-dependent Hubbard parameters obtained from band- 
structure calculations for two-color superlattices with commensurate 
(black symbols) and incommensurate wavelengths (gray symbols) 
for S2 = 10 and si = 1. Lines to guide the eye. 



of the energy scales reveals that the weak secondary laser with 
si = 52/IO considered in this example is sufficient to create 
a superlattice with e max > U. For the following calculations 
we consider a parameter range S2 = 2 to 16 and si = to 2. 
For fixed si = 0, i.e. in a homogeneous lattice, the variation 
of S2 covers a range from U ' / J ~ 1 to 60 with e max /J = 0. 
The variation of si spans a range up to e max / J ~ 4 at 52 = 2 
and up to e max / J ~ 180 at 52 = 16 with U/J being almost 
independent of si . Throughout this work the bars indicate av- 
erages over all lattice sites. 

Density Matrix Renormalization Group & Observables. 
Based on the Bose-Hubbard Hamiltonian with the site- 
dependent parameters discussed above we solve the many- 
particle problem using the density-matrix renormalization 
group (DMRG) approach lO, fl8lk which is one of the most 
powerful methods for solving one-dimensional problems of 
this type. 

The so-called infinite- size DMRG algorithm is based on an 
iterative growth of the lattice system. We start from an initial 
block of lb lattice sites with up to N\y particles per site de- 
scribed in a Fock space T\> of dimension D\>. To this block an 
additional site with up to 7V S particles described by a space T$ 
is attached to build the system space jF sys = :F b JF S of di- 
mension D sys = D\yD s . In order to mimic the thermodynamic 
limes, the system is coupled with an analogously constructed 
environment yielding the superblock JT super = jF sys <g> jF env 
with the constraint of fixed total particle number. The ground 
state |^o) is obtained by diagonalizing the superblock Hamil- 
tonian. A reduced density-matrix is formed by tracing out the 
environment, p re d = Tr env |^o)(V ; o| , and the eigenvectors 
of Pred for the largest eigenvalues are used to span the Fock 
space T\> of a new block of length 1^ = 1^ + 1 for the next 
iteration. These eigenvectors also define a non-unitary trans- 
formation matrix O which is employed to construct the new 
block Hamiltonian = O^H sys O. This cycle is repeated 
until the final length of the lattice is reached. 

In non-homogeneous systems only at the very last step of 
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the algorithm the full information about the lattice topology 
is included. Therefore, the infinite-size algorithm alone is 
not sufficient to provide an accurate description of the ground 
state for superlattices. During a second stage of the calcula- 
tion we thus resort to the so-called finite- size DMRG algo- 
rithm. Keeping the size of the superblock fixed at its final 
value, the size of the system is increased at the expense of 
the environment and vice versa. While sweeping the bound- 
ary between system and environment forth and back using the 
same operations discussed for the infinite- size algorithm, the 
Hamiltonian always takes the whole lattice into account. 

For the following calculations we use the inifinite-size algo- 
rithm for the initial phase and perform three complete sweeps 
of the finite- size algorithms in the second phase. To warrant 
convergence of the whole procedure, we increase of the block 
basis dimension D b until all observables stabilize. We typ- 
ically use D\y = 56 with up to five particles per lattice site 
TVb = N s = 5. This restriction is sufficient because we 
are most interested in the strongly correlated regime. Among 
the observables we analyze are the energy gap AE between 
ground and first excited state, the condensate fraction f c given 
by the largest eigenvalue of the one-body density matrix, the 
maximum fluctuation cr max of the occupation numbers across 
the lattice, and the visibility v of the interference fringes ex- 
tracted from the quasimomentum occupation numbers. The 
latter is of particular interest as it is directly accessible in ex- 
periment. The detailed definition of those observables is dis- 
cussed in Refs. lEltllH. 

Commensurate Lattice. We first consider the commensu- 
rate superlattice (A2 = 800 nm, Ai = 1000 nm, = 7r/4) 
with the site-dependent Hubbard-parameters depicted in Fig. 
[2] Already in Fig. [T]we have used this superlattice to illus- 
trate the difference between a generic phase diagram spanned 
by the the Hubbard parameters U/J and e max / J (neglecting 
the site-dependence of U and J) and the experiment- specific 
phase diagram spanned by the laser intensities S2 and s\. In 
both cases we show the energy gap obtained from DMRG cal- 
culations for / = TV = 30. Because the variation of S2 and s± 
affects all site-dependent Hubbard parameters simultaneously, 
the (s2j si) phase diagram is distorted relative to the generic 
(U/J, e max / J) phase diagram. Nonetheless, all relevant quan- 
tum phases can be reached by variation of the laser intensities 
alone, i.e., there is no need to vary other experimental param- 
eters like the interatomic scattering length. 

A more detailed analysis of the phase diagram as function 
of 52 and si is given in Figs. Oa) to (d), where we depict the 
energy gap, the condensate fraction, the maximum number 
fluctuations, and the visibility, resp., obtained in DMRG cal- 
culation for / = N = 30 with the commensurate superlattice. 
The superfluid (SF) phase is signaled by a small or vanishing 
energy gap, a large condensate fraction, large number fluc- 
tuations, and maximum visibility. Although the most strin- 
gent order parameter for the SF phase — the superfluid frac- 
tion lfl9ll — is not computed here, this signature allows us to 
identify the SF phase in the regime of small 52, roughly up to 
82 < 6 for all si < 2. In this regime the superlattice is shal- 
low such that tunneling dominates over on-site interactions, 
which is a prerequisite for the long-range coherence present 



in the SF phase. For si = the ratio U/J is about 1 at 52 = 2 
and reaches about 4.5 at S2 = 6. Around this value of U/J 
the phase transition from superfluid to Mott insulator is ex- 
pected in one-dimensional Bose systems in homogeneous lat- 
tices EHIiHL which is consistent with our observation. Even 
in the presence of the secondary laser, i.e. for < si < 2, 
the intensity S2 = 6 leads to U/J ~ 4.5, which explains the 
presence of the SF phase in the whole range of s\ considered 
here. 

If we increase S2 at fixed s\ = 0, the system enters the ho- 
mogeneous Mott-insulator (MI) phase, which is characterized 
by a large energy gap, a small condensate fraction, small fluc- 
tuations, and a minimal visibility. These signatures are clearly 
visible in the Fig. [3] for large values of S2 and small values of 
s\. At si = and S2 = 16 the ratio U/J is about 60 and the 
system is deep in the MI regime. 

If we now increase si at fixed S2 = 16, the modulation 
of the site-dependent Hubbard parameters grows rapidly. Al- 
ready at si « 0.6 the spread of the on-site energies be- 
comes comparable to the average interaction strength, i.e., 
Cmax/j ~ U/J ~ 60. Thus it becomes energetically fa- 
vorable to move particles from the lattice sites with largest 
on-site energies to the sites with the lowest on-site energies, 
thus creating doubly occupied sites. In this way the homoge- 
neous MI phase is broken up and the transition to a quasi Bose 
glass (BG) phase is observed. In a truely random infinite-size 
lattice, the BG phase is characterized by a vanishing energy 
gap. Intuitively, this results from the continuous distribution 
of on-site energies, which allows for the construction of ex- 
cited states by redistributions of atoms to sites with infinitesi- 
mally higher on-site energies and thus infinitesimal excitation 
energies. In contrast to a random lattice, the commensurate 
superlattice exhibits only 5 different on-site energies. This al- 
ways gives rise to finite energy gaps and extended domains 
in the phase diagram, where certain occupation-number states 
dominate the ground state — two of those domains are visible 
in Fig. Oa). Only in the transition regions between those do- 
mains the energy gaps become small. In order to approach a 
more realistic BG phase, we have to consider more complex 
lattice topologies, e.g., the incommensurate superlattice. 

Incommensurate Lattice. We repeat the above analysis for 
the incommensurate lattice (A2 = 830 nm, Ai = 1076 nm, 
= 7r/3) inspired by the experimental setup of Lye et al. 
[H. These parameters lead to a modulation of the site- 
dependent Hubbard parameters with a periodicity which does 
not correspond to an integer number of lattice sites, as seen 
in Fig. [2] Therefore, the pattern of on-site energies for 
example, is not periodic anymore as in the case of the com- 
mensurate lattice. 

The resulting phase diagrams for / = N = 30 are depicted 
in Figs. Oe) to (h). Evidently, the structure and extend of the 
SF and the MI phase are not affected by the change in the lat- 
tice topology. Only the BG phase exhibits a different behavior. 
The energy gap is clearly reduced due to the irregular charac- 
ter of the superlattice. The larger number of different on-site 
energies allows for redistributions associated with lower exci- 
tation energies and reduced energy gaps. The more irregular 
the lattice the stronger the reduction of the energy gap in the 
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FIG. 3: (color online) Phase diagrams in terms of energy gap AE/ J, condensate fraction / c , maximum number fluctuations cr max , and visibility 
v obtained from DMRG calculations for N = I = 30 with the commensurate superlattice (upper row) and the incommensurate superlattice 
(lower row). The irregularities of the contour lines are due to minimal numerical fluctuations of the DMRG results. 



quasi BG phase — eventually the true BG phase in a random 
lattice would be approached. Observables like the condensate 
fraction and the visibility are still suppressed in the BG regime 
and allow for a unique distinction from the SF phase. 

Conclusions. We have studied the phase diagram of 
bosonic atoms in one-dimensional two-color superlattices 
starting from the experimentally accessible parameters of the 
optical lattice. Following a band- structure calculation to ex- 
tract the corresponding site-dependent Hubbard parameters 
the many-body problem is solved using the DMRG approach. 
Using the intensities of the two laser fields, si and S2, as con- 
trol parameters while keeping all other experimental parame- 
ters fixed, it is possible to access all relevant quantum phases 
of the system. By following a simple path in the (51,52)- 
parameter-plane, the transition from SF to MI phase (varia- 



tion of 52 for fixed si < 0.4) and the transition from MI to 
quasi BG phase (variation of si for fixed 52 > 10) can be 
studied in detail. By comparing superlattices obtained from 
lasers with commensurate and incommensurate wavelengths 
we show that the SF and the MI phase are largely insensitive 
to the detailed lattice topology. As expected, only the energy 
gaps in the quasi BG phase change systematically when going 
to a more irregular lattice topology and approach the limit of 
a gapless BG expected in a random lattice. By measuring the 
energy gap and, e.g., the visibility of the interference fringes 
it is possible to uniquely identify the three different quantum 
phases in experiment. Thus a detailed experimental analysis 
of the rich structure of the phase diagram is directly possible, 
both from the control and the diagnostics point of view. 
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